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We propose and discuss methods for detecting quasi-molecular complexes which are expected to 
form in strongly interacting optical lattice systems. Particular emphasis is placed on the detection 
of composite fermions forming in Bose-Fermi mixtures. We argue that, as an indirect indication of 
the composite fermions and a generic consequence of strong interactions, periodic correlations must 
appear in the atom shot noise of bosonic absorption images, similar to the bosonic Mott insulator 
[S. Foiling, et al.. Nature 434, 481 (2005)]. The composites can also be detected directly and their 
quasi-momentum distribution measured. This method - an extension of the technique of noise 
correlation interferometry [E. Altman et al., Phys. Rev. A 79, 013603 (2004)] - relies on measuring 
higher order correlations between the bosonic and fermionic shot noise in the absorption images. 
However, it fails for complexes consisting of more than three atoms. 



I. INTRODUCTION 

Atomic mixtures in optical lattices represent a novel 
laboratory system for the study of ultracold matter. 
They are intrinsically clean and feature a unique amount 
of tunability. Therefore they offer outstanding opportu- 
nities to investigate fundamental many-body systems ex- 
hibiting some very rich physics. A fascinating example is 
the formation of composite states, which are induced by 
the lattice and the interaction. Such states have recently 
been detected experimentally 0, 0,0)0 ■ 

Two-component mixtures of atoms of, say, sorts A and 
B provide the simplest example of such objects which 
can be viewed as diatomic AB quasi-molecules occurring 
if the A-B interaction is strong enough. It is interesting 
to note that the formation of composites is possible not 
only for attractive Q but also for repulsive interaction 
0. In the last case the pairing occurs between an atom 
of sort, say, A and a hole B of the sort B. Moreover, 
multi-atomic com plex es of the type AB„ with n > 1 can 
form as well |^ IglllOj . Even in strongly non-equilibrium 
situations quasi-bound atomic states can exist if the re- 
pulsive energy greatly exceeds a typical single particle 
bandwidth. This mechanism was first introduced for the 
point-like defects in solid He4 and realized recently 
in the optical lattice [l^ . 

In this paper particular emphasis will be put on generic 
schemes for detecting the composites in the lattice Bose- 
Fermi mixtures since such mixtures have recently been 
studied experimentally 0, ■ We will be referring to 
the composites of the type FB„ or FB„, n — 1,2,3,... 
where F stands for a fermion and B (B) denotes a boson 
(bosonic hole). However, all the conclusions can be easily 
extended to any other type of a composite. 

The quasi-molecular complexes are induced by the lat- 
tice and consequently vanish in most cases once the lat- 
tice potential is removed. Since the standard method of 
detection in ultracold gases, namely time of flight imag- 
ing, entails the elimination of all trapping potentials, the 
detection of composites with this technique is challeng- 
ing. In the example of Bose-Fermi mixtures, the forma- 



tion of such complexes — composite fermions — would 
lead to the destruction of the Fermi surface of the orig- 
inal fermions. Therefore the disappearance of the corre- 
sponding feature should be observable in the absorption 
image. However, this disappearance would by no means 
constitute a definite proof of the formation of composites 
since a variety of phenomena, e.g. heating, could also 
lead to the same observation. 

In this letter we propose a method how such compos- 
ites can be detected directly. The method is similar to 
the Hanbury Brown and Twiss noise correlation interfer- 
ometry 15] applied to thermal sources of cold atoms 
and proposed in ref. 17] for revealing non-trivial many- 
body states in ultracold atoms. Such a method has been 
successfull y im plemented in ref. 18, 19, 20] for bosonic 
and in ref. |2lL 123] for fermionic systems. 

Noise correlation spectroscopy relies on the fact that, 
while the momentum distribution itself shows no particu- 
lar features, the atomic shot noise is correlated. In order 
to introduce the extended noise correlation method with 
a simple example, consider a model system with a sin- 
gle composite (quasi-molecule) consisting of one boson 
and one fermion. The quasi-momentum of the center of 
mass of the composite is denoted by Q. When the lattice 
is switched off, the pair breaks and the fermionic and 
bosonic wavepackets expand independently from each 
other. The fermionic or the bosonic momenta distribu- 
tions np{c\x) and ^3(^2), respectively, by themselves, 
which can be measured after expansion, do not reveal 
any low energy structure in their individual momenta 
^l^ ^.■2 - they feature almost uniform backgrounds be- 
cause binding of both atoms entails a wide range of their 
relative momenta. However, momentum conservation re- 
quires the momenta to add up to Q. Therefore, if in an 
individual measurement, the fermion is found to have 
qi, the boson must have q2 = Q — qi and the correla- 
tion between the two enables one to measure the quasi- 
momentum of the composite. 

In the more general case with more than one pair in 
the system, the correlations between the bosonic and 
fermionic momentum distributions, which allow one to 
determine the quasi- momenta of the pairs, will only be 
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observable in the noise 5nB,F = ns.p — {nB,F) of the 
measured momentum distributions. This is due to the 
fact that, besides the correlated pair events, there are 
also events due to the components originating from dif- 
ferent pairs. The key point is that the momentum distri- 
bution of the (centers of mass of the) FB composites is 
reproduced by the noise correlator in the far field of free 
expansion 

Icf{Q) = J dX2{SnF{q,)SnB{q2)), (1) 

with Q = Qi + q2 fixed, and qi = rapl^i/t = Q — q2, 
q2 = 1713^2 /t being the dependent and independent 
variables, respectively. Here mp, ttib are the masses of 
fermionic and bosonic atoms, respectively, and Xi (X2) 
denotes the position in the far field where a fermion (bo- 
son) is detected. We point out that ^ is essentially 
the Fermi distribution function pcf(Q) of the composite 
fermions ( pcpiQ) w 1 for |Q| < Qcf and Pcf(Q) ~ 
for IQI > Qcf, where Qcf stands for the Fermi surface 
wavevector of the composites) . We also show that in the 
repulsive Bose-Fermi mixture, the role of the center of 
mass momentum Q is played by Q_ = qi — q2 similar 
to the case of the repulsive two-component bosonic mix- 
ture These properties of the correlator will be 
discussed in detail in the appendix 1X1 and extended to a 
general case of FB„ and FB„ complexes in section III CI 
Detection of the Fermi edge in IcpiQ) may require 
quite low temperatures T determined by the Fermi en- 
ergy EcF of the composites. However, since a typical 
binding energy Ei, of the composites can be much larger 
than EcFi a signature of the composites will be seen 
as long as T < Eh (for T > Ejj, the correlator ^ is 
essentially zero). Indeed, the very fact that Icf{Q) is fi- 
nite and proportional to a typical density of the original 
fermions or bosons is a direct indication of the presence 
of the composites. Accordingly, as T is lowered below 
EcF, the formation of the Fermi edge in Icf{Q) will be 
observed. 

In principle, detection of multi-atomic complexes — 
quasi-molecules with the number of the atomic con- 
stituents exceeding two — can be done by extending the 
noise-correlator method, as explained in Sec lll CI How- 
ever, as will be discussed in Sec lII Dl such an extension to 
higher order correlators can lead to a dramatic increase 
of uncorrelated noise. For this case, we propose an indi- 
rect method to observe the complexes (Sec lIII|l . In this 
method, the presence of complexes of the type FB„ is re- 
vealed by imaging the density-density noise correlations 
of the original bosons, if their number is incommensurate 
with the number of the lattice sites. Then the resulting 
image is similar to the one observed in a bosonic Mott 
insulator 0, 0| . Such a feature (at zero temperature) 
turns out to be a consequence of a strongly interacting 
many-body ground state in which individual momenta of 
bosons are poorly defined and, therefore, strongly fluctu- 
ate. 



II. DESCRIPTION AND DIRECT DETECTION 
OF COMPOSITES IN OPTICAL LATTICES 

A. The nature of the composite pairs 

Here we will give a brief account of the nature of 
equilibrium quasi-molecular states. In general, quasi- 
molecular (composite) states form when the binding en- 
ergy Eb exceeds the temperature T as well as the gain in 
kinetic energy the constituents would encounter by de- 
localizing from each other. For attractive interactions 
between bosons and fermions, composites states of the 
type FB are formed. For repulsive interaction, on the 
other hand, the quasi-molecules represent bound states 
of one fermion F and a bosonic hole B or, similarly, one 
fermionic hole F and a boson B. Binding of one fermion 
with more than one boson or bosonic hole is also possible 
if the number of bosons per site is large enough and the 
magnitude of the Bose-Fermi interaction |?7Bi?| substan- 
tially exceeds the repulsion between bosons Ub b [io| . In 
this event the composites FB„ (for Ub _f < 0) or FB„ (for 
Ubf > 0) with the integer n « \Ubf\/Ubb can form. In 
recent experiments .ISi. il4j complexes of the type FB„ 
with n as large as n = 3 may have played an important 
role. 

Depending on the parameters, the effective interaction 
between composite fermions can be tuned to become at- 
tractive or repulsive. In the case of effective attractive 
interaction, a p-wave superfluid of pairs of composite 
fermions will form For effective repulsive inter- 

actions and low filling of composite fermions the ground 
state is a Fermi liquid while near half filling insulating 
states may appear [Toj |. For clarity, we will limit our- 
selves here to the case of composite fermions forming a 
weakly interacting Fermi liquid with a well defined Fermi 
surface. 

A pure quasi-molecular phase FB„ will only exist when 
the number of bosons N b and fermions Np is commen- 
surate, i.e. Nb — nNp. Conversely, in the case of a 
pure phase of composites of the type FB„ the number 
of bosons must be Nb — N ■ k — Np ■ n, with the inte- 
ger k > n [23 and N standing for the total number of 
the optical lattice sites. Free fermions or bosons will be 
present when these conditions are not fulfilled, or when 
going away from the strongly interacting regime, where 
the composites can partly dissociate. Then two types of 
Fermi surfaces can coexist — one for the composites and 
the other for the original fermions. Moreover, a portion 
of the B-bosons can form a BEC l23| . 

Here we will consider the limit of low density of the 
composites, so that the size of one molecule am is much 
less than the average distance between them. In this 
case the description becomes quite straightforward in 
terms of the centers of mass of the composites only. As 
o~m/fo increases and eventually becomes comparable to 
1 , the composite fermions will dissociate according to the 
scenario l24l. 
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B. Detection in weakly interacting Bose- Fermi 
mixtures 

In the weakly interacting Bose-Fermi mixture with- 
out composites present, practically all important infor- 
mation can be deduced by imaging the bosonic nB(q) 
and fermionic np(q) quasimomentum distributions. In 
this method, the quasimomentum distribution is first 
mapped to the momentum distribution by removing the 
lattice potential on a time scale which is slow with 
respect to single particle physics but fast compared 
to quasimomentum changing collisions. After subse- 
quent free expansion for a time t the quasimomen- 
tum distribution is detected by imaging the density 
|25ll26|. Apart from unimportant multiplicative factors, 
nsiq) « / dxexp(-zqx)(V'|j(x)V'B(0)) and npiq) w 

J dxexTp{—iqx.){il;p{-x)ipF{0)), where V'-b(x) and V'f(x) 
are the bosonic and fermionic field operators, respec- 
tively, just before the lattice was removed and the av- 
eraging (...) is taken over the corresponding state in the 
lattice. In the weakly interacting regime, the bosonic 
quasi-momentum distribution reflects the off-diagonal 
long range order (ODLRO) inherent in the one-particle 
density matrix pb(x, x') = (?/)^(x)f/)B(x')), resulting in a 
very small spread of quasi-momenta. The fermionic one- 
particle density matrix pi?(x, x') — {tpp{x.)tpFi^')) fea- 
tures Friedel oscillations which are due to the sharp edge 
in the Fermi distribution np{q) (at T = 0): npici) ~ 1 
for |q| < qp and np{q) « for |q| > qp, where qp stands 
for the Fermi wavevector. 



C. Direct detection of composites 

The method described above is not applicable to com- 
posite fermions, since they typically dissociate during the 
expansion. However, the center of mass momentum Q of 
a FB„ atomic complex is conserved during time of flight. 
Therefore the momenta of all its n-t-1 constituents after 
expansion are directly correlated and add up to Q. We 
point out that the observation of such correlations consti- 
tutes direct evidence for the existence of the composites. 
In the following we will describe how these correlations 
can be used to directly detect the (quasi)momentum dis- 
tribution of the composites in the lattice. 

For an ultracold gas that contains many composites, 
the fermionic (bosonic) momentum distributions, i.e. the 
absorption images after time of flight, represent the mo- 
menta of many fermions (bosons). Consequently, the cor- 
relations between the constituents originating from one 
composite are not easily visible, since it is not clear which 
of the many bosons and fermions made up one composite. 
The correlations can therefore not be determined from 
the form of the momentum distributions but only from 
the atomic shot noise visible in them. Since the correla- 
tions are between one fermion and n bosons, the corre- 
sponding correlator is {dnp{'X.i)SnB{'^2)---SnB(^n+i)) ■ 



To determine the amount of correlation, this correlator 
must be integrated over n + 1 coordinates under the con- 
straint that the center of mass momentum Q of the con- 
stituents is fixed: 

-?^n+i(Q) = y"dq2...rfq«+i (2) 

(: Snp{c[i)SnB{q2)---SnB{<ln+i) ■), 

with Q = qi + q2 + ... + q„+i and qi = mplLi/t, q2 = 
7715X2/^, q„+i — ruB^n+i/t. Here : ... : implies the 
normal order of the operators — when all the creation 
operators stay to the left of the annihilation ones. The 
averaging (...) is performed over a particular realization 
of a many body eigenstate j^*) and, then, over many real- 
izations. In appendices A and B we show that the amount 
of correlation /„+i(Q) is proportional to the number of 
composites with given quasi-momcntum Q, that is 

In+i{Q) - Pcf{Q) - C, (3) 

where pcf(Q) — (c^(Q)c(Q)) and c'''(Q),c(Q) are the 
composite fermion creation and annihilation operators 
in the momentum space and C is a constant term. To 
summarize, the quasimomentum distribution pcf(Q) of 
the composites can be found by measuring the correla- 
tions between experimental absorption images according 
to the equation (jSJ. 

The relation is proved rigorously in appendices A 
and B. It is instructive to follow the main steps of the 
proof: The operators 

Snp{X) = ^/jl{X)i;p{X)-{^\^l{X)^p{X)\^),{i) 
SriBiX) = Vt^(X)V^B(X)-(vl/|^t^(X)^B(X)|M'),(5) 

appearing in equation (|2Jl are expressed as functions of 
the onsite fermionic and bosonic creation and annihi- 
lation operators with i being the site in- 
dex, in the lattice before expansion with the help of 
eas. (jA4IA5IA3|) . Inserting these expressions into eq.Q 
yields sums over terms with 2{n + 1) creation and anni- 
hilation operators on different sites. Many of these can 
be set to zero, resulting in 

/„+i(Q) ^ 5^e*'^(''--^'([clc,-%clc.]), (6) 

where c|, Ci are the onsite composite fermion creation- 
annihilation operators as defined in ea. (|Bl(l of the Ap- 
pendix ^ and Xi denotes the site i coordinates. The 
first term in brackets in eq.® is pcpiQ)'- as can be no- 
ticed, after introducing the composite Fermi operators 
in the Fourier space c(Q) ~ exp(— jQxi)ci, one finds 
equation (PJ. While pcf{Q) gives the occupation num- 
ber of composites with quasimomentum Q, the second 
term is the constant C which takes care of the normal- 
ization J (iQ/„+i(Q) — 0. The last relation follows from 
the observation that the integral J dQ /„+i(Q) is exactly 
{6Np : (SNb)" '■) where SNb.p denote fluctuations of the 
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total numbers of bosons and fermions Nb,f, respectively. 
Since Np is diagonal on any many-body eigenstate |^), 
such a mean is exactly zero for any particular realiza- 
tion of l^f), that is 5Np\'^) — 0. Both contributions in 
eas. H3l6|l can easily be distinguished because they have 
very different structures: the first one is concentrated 
within the Fermi surface of the composite fermions and 
the second is spread uniformly over all momenta in the 
first Brillouin zone. 

In the case of the complexes FB„ ( see the discussion 
in Appendix IXlbelow eas. l|A9IA10|l ') the only change is in 
the definition of the momentum Q in eqs.Q. It becomes 
Q = qi - q2 - . . . - qn+Lj 

As discussed in ref.[2J|, decreasing the F-B interac- 
tion Ubf will eventually result in coexistence of three 
types of quasi-particles: some portion of free B-bosons 
and F-fermions and the composites. In order to de- 
tect the composites in this case, the noise must be mea- 
sured on top of signals from the structured backgrounds 
(ni?(q)), ('^^(q)) in the fermionic and bosonic channels. 



D. Correlated versus uncorrelated noise 

Here we discuss the detection requirements under the 
assumption that the atomic quantum and shot-to-shot 
noise is dominant. The above results indicate that mea- 
suring the noise in the density correlator of the n+1 order 
can provide a one-particle density matrix for the centers 
of mass of the composites. It is important, however, that 
the uncorrelated noise due to photon shot noise, inac- 
curacy A in determination of the number of particles in 
each detection bin, etc., is much less than the correlated 
noise in one shot. Otherwise, many shots S are required 
in order to reduce the effect of the uncorrelated noise (by 
the factor ~ 1/^/S). Here we will show that, while for di- 
atomic composites the uncorrelated noise (with respect 
to the correlated one) is insignificant for large enough 
numbers Np ~ Nb (decays as l/y/Np), for multi-atomic 



with n = 3, .. it grows as 



N 



n/2-l 



composites FB^ 
with the case n = 2 being marginal. This makes the di- 
rect scheme of measuring the density matrix of the n > 2 
composites impractical because the required number of 
shots grows as iS* ~ Np~'^. This effect is caused by a 
strong increase of the phase space with n when the com- 
posites are dissociated as a result of the lattice release. 

In what follows we will normalize the correlated and 
uncorrelated noises by the corresponding measured signal 

Isig= J dX2...rfX„+i(nF(Q - q2 - ••■ - qn+i)) 

("s(q2)>--("s(qn+i)) (7) 

where q/c — niB^k/t. This integral can be estimated as 
Isig ~ N^Np/Nhin = "-"^F^V^bin, whcrc wc considcr 
exact matching Nb = nNp and A^bin stands for the total 
number of the detection bins. In the following estimates 
we will be using the relations Nb ~ Np ~ N, so that the 
actual numbers of particles can be replaced by the total 
number of sites (provided the densities are finite). 



The uncorrelated noise is caused by technical inaccu- 
racies A(Xi) in the detection of the number of, e.g., 
fermions in each detection bin. The effect of such noise 
becomes critically magnified for large n even under the 
condition that |A(Xi)| is much smaller than a typical 
quantum fluctuation ^ y/ Np / Abin ^ y/ A/A"bin S> 1 of 
the number of detected particles per each bin. We note 
that, under the ideal condition A(Xi) — 0, quantum fluc- 
tuations of the particles number in each bin do not give 
rise to any uncorrelated noise in the correlator This 
is due to the fact that the center of mass momenta of the 
composites are good quantum numbers - in contrast to 
the momenta of the constituents - and therefore the cor- 
responding quantity, the density of the composites rep- 
resented by the operator in eq.ld)) does not fluctuate in 
any particular many-body eigenstate \'^). If, however, 
the measurement is done with an error A(Xi), the opera- 
tor determining the correlator Q acquires the contribu- 
tion AI = J dX2...(iX„+iA(Xi)(5ns(X2)...(5nB(X„+i), 
which can be considered as a fair estimate of the 
total error under the condition |A| <C \/Np/Nhm- 
Clearly the mean of it is zero, (A/) — 0. How- 



ever, its rms fluctuation E = \J {AIAI) / Igig (normal- 
ized by Isig) is not zero. We estimate: (A/A/) ^ 
(A2)(/dXi/dX2(<5nB(Xi)nB(X2)))" ^ (A2)A", 
where we have taken into account that the quantum 
fluctuation of the number of bosons in each bin is ~ 
N/Nbin and we have considered no correlations be- 
tween the bins (and omitted factors dependent on n like 
n!, n", etc.) . Given this estimate and the one for I^ig in 
eq.lUI), we find 



E : 



(|A|)iVbin 
Nl+n/2 



(8) 



where (|A|) stands for the rms- fluctuations of A(Xi) . 

The quantity E charaterizes the uncorrelated noise 
contribution in one shot. If S shots are performed, it 
is reduced as i? — > E/^/S so that, in principle, the 
effect of the uncorrelated noise can be reduced below 
the correlated contribution In+i/Isig after taking a suf- 
flcient amount of the shots. In reality, the required S 
must scale not faster than S ~ C(l) with Np. The 
correlated contribution ((BJ is given just by the total 
number of the composites ~ Np in each bin: In+i ^ 
Np/Nun ^ N/Nbin. Hence, In+i/Lig ~ iV"" and 
A///„+i - (|A|)A^6,;„iV"/2-7Vs^. Thus, the required 
number of shots is 



|ApA4 



(9) 



For n = 1 , this condition can be achieved for large enough 
A^ and not very large number of bins A'bin- For n = 2, 
the condition is marginally reachable, and, for n > 2, it 
becomes impossible to fulflll. 

It is worth mentioning that, in the case of composite 
bosons, say, AB„, where A labels some boson of a kind 
different from B, the condition © becomes less strict. In- 
deed, 1,1+1 in this case acquires an additional factor 
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~ N due to the ratio of the correlation volumes of a bo- 
son and a fermion. Indeed, the sum J^ij e**^^^'"''^^ (c|cj) 
can be estimated as ~ N£^^, where ^ is a typical dis- 
tance on which the correlator (c|cj) becomes zero. For 
fermions it is given by the interparticle separation Tq. For 
bosons in BEC, it is the size of the lattice and one finds 
In+i ^ N'^/Nbin- Accordingly, keeping the signal and 
the uncorrelated noise estimates the same, we find the 
number of the required shots S » \A\'^N^-^N''-'^ . This 
implies that for composite bosons the method described 
above can be used if n < 4, that is the total number n+1 
of the constituents should not exceed five. 

Concluding this section, we mention that performing 
shot-to-shot averaging should be done with care, be- 
cause, otherwise, the noise scaling with TV can become 
even stronger. The operators of fluctuations 5np — 
np — (np), 5nB — ub — {nu), ^iia — ua — {tia) are de- 
fined with respect to averaging over a particular ground 
state characterized by given total numbers of atoms 
Np, Nb, Na- In other words, the average number of 
atoms (np), {tib), («a) per each bin must be determined 
in each shot. Then, the final value In+i can be aver- 
aged over the shots. Otherwise, shot-to-shot fluctuations 
of Np^B,A which typically scale at least as ^JNp^B,A^ 
respectively, will introduce stronger noise if the means 
{np), (tib), {tia) are understood as the total ensemble 
means. In the latter case the difference between the mean 
in one particular realization and the total one serves as 
an additional uncorrelated noise. Then, one can esti- 
mate A ~ VN/Nf,in and, accordingly, find A///„+i ~ 
jV(n-i)/2^y^^ so that the required number of shots scales 
as ~ N'^~-^, implying that the detection of the com- 
posite fermions with n > 1 is not possible. Similarly, for 
the composite bosons: A///„+i ~ N^"-^^'>/'^/^/S, and 
the detection for n > 3 becomes impossible. 



III. INDIRECT INDICATION OF THE 
COMPOSITE FERMIONS IN THE BOSONIC 
DENSITY-DENSITY NOISE CORRELATOR 

Here we discuss the situation when all bosons are 
bound to be the constituents of the composite fermions 
so that there is no ODLRO in the lattice. The question 
is how the bosonic density-density correlator 

hbiY) = y"dX(*| :nB(X + Y/2)nB(X-Y/2) : I*), 

(10) 

(taken in the normal order) is affected by such state. 
Below we show that (|10() exhibits the Bragg structure 
similar to that observed in the Bose Mott insulator (MI) 
[Tsl Hfll l . The relations (|B6IB7(I in the physical space of 
the composites (IB2p allow expressing the correlator (jlOf) 



as the mean over the operator 

/,,(Y) = n2^(H-e^(™«Y/*)(-'-''^))clc]c,c, 

+ n(n-l)^clc„ (11) 

i 

where the last term, which is ~ o{l/N) with respect to 
the first one, can be ignored. Taking the means and re- 
placing (clcjCjCi) = v'^ {1 — 6ij) , with 1/ being the average 
onsite population of fermions, we find that ea. lll|) has ex- 
actly the same structure observed for the bosonic Mott 
insulator the term ~ expjt (m f. b Y /f ) (x^, — xj)] 

is the correlated contribution 'l5'| which peaks as long 
as Y matches the positions of the Bragg peaks (in other 
words, Y = modulo the primitive vectors of the recip- 
rocal lattice). It is important that this structure exists 
for any filling factor of the composites i' < 1, that is, it is 
not a sole property of the Mott insulating ground states. 
Regarding the fermionic correlator 

Iff{Y) = J dX(*| : np{X + Y/2)np{X - Y/2) : |*), 

(12) 

we note that it must exhibit the antibunching peaks, 
that is, ///(Y) — for Y = {modulo the primitive 
vectors of the reciprocal lattice) regardless of the na- 
ture of the ground state because the normal product 
: np{'X.)np(X.) : is zero (due to the Pauli principle and 
the periodicity of the lattice). This can be verified di- 
rectly with the help of eas. HA4IA3l) . Then, Eq.lO be- 
comes a periodic function J//(Y), with the periods de- 
termined by the primitive vectors of the reciprocal lat- 
tice, and one finds ///(O) ^ / dc[J2ijki exp(iq(xi -I- x/j — 

Xj — xi)){fjflfifj) = due to the Fermi statistics: 
fifj+fjfi^O. 

The physical reason why the bosonic correlator H1U|) 
shows bunching peaks (|ll|l in a strongly interacting mix- 
ture stems from a very general principle, which is best 
explained by starting from the opposite case. In any 
non-interacting pure quantum mechanical state the cor- 
relator (\['o|(5riB(X)(5riB(Y)|^'o) (not taken in the nor- 
mal order!) must be exactly zero. To understand 
this one must remember that the density operators in 
the far zone are equivalent to the (quasi-)momentum 
operators in the lattice and that non-interacting sys- 
tems are characterized by well defined quasi-momenta. 
Thus, there are no fluctuations in each particular many- 
body (non-interacting) eigenstate j^I'o) and the correla- 
tor (5'o|(5ns(X)(5ns(Y)|^'o) is zero. In other words the 
operator is zero on space of states |5'o) of the ideal 
systems: that is, (5nB(X)|5'o) = 0. In the thermodynam- 
ical limit at low T the normal ordering in ea. HlOl) can 
be ignored for bosons (as opposed to fermions!). Conse- 
quently, a weakly interacting BEC of B-atoms will lead 
to essentially uniform distributions in ea. HlU|) |2^. 
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Meanwhile in a quasi-molecular phase of strongly in- 
teracting Bose-Fermi mixtures the momenta of the orig- 
inal bosons are undefined in each pure many-body state 
I^E") in complete analogy with the MI. Accordingly, the 
densities in the far zone fluctuate strongly and produce 
the bunching behaviour in (|If |l . In this case quantum 
fluctuations are responsible for the Hanbury Brown and 
Twiss effect. Setting it more technically, in each many- 
body eigenstate the correlator (\E'|(5n_B(X)5n_B(Y)|\E') is 
proportional to the interaction vertex and, if the 
interaction is small, the correlator is, practically, zero. 
In general, having a strong interaction regime is the re- 
quirement for observing the distinct patterns in pure 
states. This condition is more general than the one for 
having energy gap with respect to single-particle excita- 
tions 18]. In the state mentioned in Sec III CI where all 
three types of quasi-particles are present, such a gap is 
zero and the composites can still be well defined. Accord- 
ingly, the Hanbury Brown and Twiss pattern should 
be observable. Obviously, some care should be taken in 
order not to confuse the pattern due to the composites 
with that caused by thermal fluctuations 

IV. CONCLUSION AND 
ACKNOWLEDGEMENT 

We have suggested the extended density-density corre- 
lator method for the direct imaging of the composites in 
optical lattice and have evaluated its feasibility for multi- 
atomic fermionic and bosonic complexes. The method is 
expected to work well for fermionic composites consisting 
of less than four atoms and for bosonic ones consisting of 
less than six atoms. Above these numbers, the uncorre- 
lated noise becomes too large. 

As an indirect method for detecting the composites, 
measuring the bosonic density-density correlators must 
reveal the typical Hanbury Brown and Twiss structure 
first observed for the bosonic Mott insulator. We point 
out that this observation is a generic consequence of 
strong interactions in a two-component mixture in op- 
tical lattices. 
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APPENDIX A: STRONGLY BOUND FB AND FB 
PAIRS 

Here we will give a more extended treatment of the 
FB case discussed in the introduction. It is important 



that the correlator |^ turns out to be proportional to 
the Fourier representation of the in-situ two-body den- 
sity matrix (V'l^(xi)i/'B(x2)V'F(yi)'0B(y2)) with respect 
to the coordinate Xi — yi. As a matter of fact, since a 
boson can only be found close to a fermion, one can set 
xi = X2, yi = y2 in this correlator without affecting its 
low energy properties (as long as am/fo <^ 1). 

In the following we will justify the statement /c-f(Q) ~ 
Pcf{Q)- We assume that the wavefunction of relative 
motion $(r) of a boson and a fermion in a pair is of the 

s-wave type $(r) « exp(— |r|/am)/am^ with am being 
of the order of the lattice constant. [In the case of the 
non-s-wave BF pairing, the correlator vanishes after 
integration over, say, q2 with Q kept constant. We do 
not consider such an exotic possibility]. We will see in 
the following that once IQa^] <^ 1, the extent of the 
wavefunction a„i has no effect on /c_f(Q) or pcf(Q)- In 
order to calculate /cf(Q), we expand the fermionic and 
bosonic operators in the optical lattice as 

i^f(x) = ^I4^o(x-x,)/., (Al) 

i 

V'b(x) = ^Wo(x-x,)6., (A2) 

i 

where Wo{x—Xi) stands for the Wannier function located 
at ith site. For simplicity, we consider VFo(x) to be the 
same for fermions and bosons. Upon free expansion of a 
particle of mass m, the Wannier function in the far zone 
becomes (apart from a numerical coefficient) 

iyt(X-x,:) = (^!!^^'^%-»qx.+™x=/2t yj/(q) 

(A3) 

where q = m~K./t (we employ units in which h — 1) and 
PF(q) = J dx' exp(— iqx')VFo(x') stands for the Fourier 
transform of the Wannier function. In the far zone the 
operators I|A1IA2|I become 

V'f(X) = ^m(X-xO/,, (A4) 

i 

V'b(X) = Y,WtO^-^^)h, (A5) 

i 

where the corresponding mass mp, Tas is replacing m. 

Then, substitution of eas. ljA4IA5|l into the definition 
of the densities in the far zone (|4I5|I and employing them 
in eq.QJ give 



7 



IcFiQ) = / d(l2 



[(/^(Q-q2)&^(q2)/(Q-q2)6(q2)) 



(/^(Q-q2)/(Q-q2))(6t(q2)6(q2))] 

I 



(A6) 



where /(q) = ^ exp("iqx,)/, and 5(q) = 

Ej exp(— iqxi)6i are the Fourier representations of 

the onsite operators. We note that /cf(Q) scales with 
the expansion radius R ^ t as ~ that is, as though 

it is just a density of some particle undergoing free ex- 
pansion. 

The momenta involved in W^(q) are of the order of the 
inverse lattice constant ai. In fact, as the lattice is being 
ramped down adiabatically with respect to the single par- 
ticle states, the Wannier function undergoes a transfor- 
mation from that corresponding to the deep lattice to the 
one in the very shallow lattice. In the latter case, VWx) 
can be found explicitly from the definition (1) of ref . |29l| . 
where for the Bloch function one can use just the expo- 
nential exp(ikx). Then, the square of the Fourier trans- 
form |M^(q)p trivially becomes the volume of the elemen- 
tary cell ^Ibz ~ o,f when q belongs to the first Brillouin 
zone and zero otherwise. Accordingly, the integration 
/ dq2... in ea. HA6|l proceeds over the first Brillouin zone 
/g^dq2... with the factor |l^(Q - q2)W^(q2)P = 1 as 
long as the typical values of Q are much smaller than the 
inverse lattice constant ai (for Q being outside the first 
zone, this factor is essentially zero, if the Fermi "sphere" 
has a volume w Qcf much less than l/flsz)- In par- 
ticular, this implies that the "imaged" composites are 
localized mostly within the first Brillouin zone similarly 
to the rampdown procedures [2^ used for the atomic 
imaging. As can be seen, the spreading over the next 
zone is suppressed by the factor Q^pilsz <C 1 due to 
the small fraction of the phase volume in the integral 
(|A6|) which allows access to the center of the next zone. 

If bosons and fermions do not bind together, the corre- 
lator in ea. (jA6|) is essentially zero. If, however, there are 
tightly bound composites FB in the lattice, the correla- 
tor become finite whenever a boson is found within the 
proximity am of a fermion. This is indicated in a choice 
of the physical space 



ii,{"Q}>^n(^Q)"''io)' 



(A7) 



Q 



= ydq$(q)/t(Q/2 + q)6t(Q/2-q), (A8) 

with riq = 0, 1, of weakly interacting composite fermions 
(in the considered limit Vo ^ am). It is represented by 
the composite creation-annihilation operators Cq, cq of a 
BF pair with total momentum Q in terms of the Fourier 
components of the fermionic and bosonic operators and 
of the relative wavefunction $(q). It is important that 
the low energy properties do not depend on details of 



$(q) in the limit of low density of the composites. This 
can be seen by direct calculation of eq. ljA6|l over the state 
(|X7|) . Specifically, the state /(Q - q2)6(q2)|l, {nq}) 
contains contribution $(Q/2 — q2)|0) coming from the 
action of /(Q — q2)6(q2) on the term describing just 
single composite Cq|0) without affecting the others. It 
also contains the exchange terms involving pairs of the 
composites with different momenta. These exchange 
terms, however, can clearly be neglected because they 
contribute factors {am/ro)^ <C 1. Thus, neglecting these 
terms, we find /g^ dq2 ({"-q}, 1|/^(Q - q2)&^(q2)/(Q - 
q2)6(q2)|l,{nQ}) = /3^dq2|$(q2)P riQ 
({nQ},l|c^(Q)c(Q)|l,{nQ}), Since /^^ dq2|l>(q2)P = 1 
as the normalization condition. Thus, in the phys- 
ical space (|A7IA8|I the operator Jg^ dq2/t(Q — 

q2)6'''(q2)/(Q - q2)&(q2) is equal to c'I^(Q)c(Q) and 
its averaging produces the momentum distribution 
PCFiQ) = ({nQ},l|ct(Q)c(Q)|l,{nQ}) of the com- 
posite fermions as discussed in the Introduction below 
eq.©. 

As can be directly checked, the uncorrelated contri- 
bution /^^dq2(/t(Q - q2)/(Q ~ q2)) (6^ (q2)fe(q2)) in 
ea. (|A6|) is independent of Q as long as |Q| ^ and, 
thus, can be replaced by a constant C. This constant 
can be restored from the normalization on the density 
of fermions (and bosons). Finally, the expression be- 
comes /c_f(Q) ~ Pcf(Q) — C, which is a particular case 
of a more general eq.® discussed in the section HTni 

We also comment on the case of purely repulsive inter- 
actions when the pairing occurs between fermions and 
bosonic holes 0- The corresponding physical space, 
then, becomes 



\h{nQ})=Yl{clr'^Y[ibir\0) (A9) 

Q i 

= ydql>(q)/t(Q/2 + q)6(q-Q/2), (AlO) 



where fc = 1, 2, .. and $(q) is now a wavefunction of rela- 
tive motion of a fermion and bosonic hole. Accordingly, 
the total momentum Q is now the difference of the mo- 
menta carried by the fermion Q/2-f q and by the bosonic 
hole q— Q/2. Then, as discussed in the Introduction, the 
momentum Q in eq. Q is to be replaced by the difference 
of the momenta: Q_ = qi — q2. 
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APPENDIX B: THE PHYSICAL HILBERT SPACE 
FOR STRONGLY BOUND COMPLEXES FB„ 

Insensitivity of the low energy physics to the internal 
structure of the composites discussed in the Appendix IXI 
for the case FB (as an example) allows a simplification of 
the analysis by projecting the Hamiltonian as well as any 
measurable quantity into states having either no particles 
on a site or l+n particles — one fermion and n bosons. In 
other words, the constituents of a composite are located 
at the same site and the probability to find an unbound 
boson or fermion on another site is vanishingly small. 
This means that the physical operators are 



Then, 
basis 



(Bl) 



the physical Hilbert space is represented by the 



K{M>-n(4no)' 



(B2) 



where the product is taken over all lattice sites and 
can take only two values = 0, 1. As can be easily 
checked, the annihilation and creation operators of the 
composite fermions obey the standard anti-commutation 
relation [ci,cj]+ — Sij in the space HB2II . 

The projection of the measured quantities into the 
space (|B2|) can be done by considering processes leading 
to jumps of the composite fermions only. If the compos- 
ites are strongly bound and there are no free bosons or 
fermions, this is a very accurate approximation. Then, 
we ignore blbj and fjfj acting on |n, {fii}) unless i — j 



which gives onsite local populations determined by pres- 
ence or absence of a composite on the site i. So, in general 



i"3 



nc-CjS. 



fjfj — c\ci5i. 



(B3) 
(B4) 



Concerning the quartic operators fifjb\.bi in eq.(Q 
which are important for the detection of the FB quasi- 
molecules, they must be set to zero (if acting on |1, {ui})) 
unless i = i, k = I OT i = k, j = I. For n = 1, the follow- 
ing relation holds 

fifjbibi = SikSjiil - 5ij)clcj + 5,j5kic\c^clck. (B5) 

Similar relation can be established for any value of n for 
the operator linear in /J/i and of the order n in 6^6^. We 
will not, however, present it explicitly here. Another set 
of useful relations is: 



fiflfifj = [SijSki- Sii5kj]clclckC^, (B6) 
biblbibj = n'^[Sij6ki + SiiSkj]clclckCi 

+ n{n - l)S.aSijSikclci. (B7) 



As long as the low energy properties are concerned in 
the low density limit, the introduced truncation of the 
phase-space is accurate. Similar rules can be formulated 
for the case of the FB„ complexes. In the case of n > 1 
the derivation of eq.® proceeds as follows. After sub- 
stituting expressions ljA4IA5|l into eas. H4l5|l and, then, 
into eq.Q, we arrive at the multiple lattice sum of the 
correlators in eq.Q: 



/„+l(Q) ~ J dq2...dq„+l e*[(Q-q2...-q. + l)x,i,i+q2X.,,,+...+q„ + iX,„^,,„^J 

milfn - mlfnl"^)) ■■ (bib,, - - V+i I*)) : I*) (B8) 

I 



It is enough to use the normal ordering : ... : for 
bosonic operators only because just two fermionic 
operators are involved. The terms with ii ^ ji 
which project into the one-particle density matrix 
of the composite fermions contain 2n bosonic opera- 
tors together with the fermionic pair under (^|...|^'). 
These are (*|/^/,, : blbj,...b] : I*) 



■bl,,b,^+. 

1 + 1^3132- 



where 



the 



ml{blrf3dbnr\m.^2■■■Sn^ 

'^K*|cliCjil*><5ni2--*<ii„+i'5jij2--'5jij„+i: 
definition (|B1|I is utilized. This generates the first 
term in the brackets of eq.®. The terms with 
ii — ji, after being summed, produce zero due to the 
conservation of the total number of fermions. This 



explains the presence of the last term in the brackets 
of eq.®. Other terms carry less than 2n bosonic 
operators averaged together with the fermionic ones. 
For example, one term which has 2{n — 1) bosonic 

operators is S.,,A^\blb.,\^)mflfn " (*l4/.J*» : 
bl,^bj^...bl ^jfej„_|_i : I ^), where the factor (5^2 j2 comes from 
the relation ljB3|) . It does not contribute to the terms 
with ii ^ ji because the integration over q2 in eas. lj2IB8|) 
has a structure J dq2 exp(i(Q — q2):x.i-^j-^ + ''■<i2^i2j2)^i2j2 
and, thus, selects only terms with ii — ji. The same 
logic applies to other terms with fewer number of bosonic 
operators. 
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